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Normal state spin dynamics of the recently discovered iron-pnictide superconductors is dis- 
cussed by calculating spin structure factor S(q,ui) in an itinerant five-band model within RPA 
approximation. Due to the characteristic Fermi surface structure of iron-pnictide, column like 
response is found at (n, 0) in extended Brillouin zone in the undoped case, which is consistent 
with the recent neutron scattering experiment. This indicates that the localized spin model 
is not necessary to explain the spin dynamics of this system. Furthermore, we show that the 
temperature dependence of inelastic neutron scattering intensity can be well reproduced in the 
itinerant model. We also study NMR \jT\T in the same footing calculation and show that the 
itinerant model can capture the magnetic property of iron-pnictide superconductors. 
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Recently discovered 1 ) iron-pnictide (or transition- 
metal pnictide) superconductors have attracted much 
attention and have been studied intensively since they 
show superconductivity at rather high temperatures up 
to 55K 2,3 ) at present. Novelty of this series of compounds 
is not limited to their high Tc: entangled band struc- 
ture, coexistence of hole and electron Fermi surfaces, 4 ' 5 -' 
origin of the magnetism, 6,7 ) interplay between mag- 
netism and superconductivity, 8 ' 9 ^ role of lattice struc- 
ture, 10, nodal or nodeless behavior in the supercon- 
ducting phase, 12 ~ 15 ) robustness against impurities in the 
conducting layer 16-18 ) and so on. All of these features en- 
rich the physics of this compound. In this paper, we pay 
attention to magnetic properties of this series of com- 
pounds, especially the normal state spin dynamics. 

For the study of spin dynamics, NMR and inelas- 
tic neutron scattering experiments are powerful tools. 
There already exist systematic works of NMR, 12,19 ) and 
some results of inelastic neutron scattering on iron- 
pnictides. 20 ~ 24 ) The normal state behavior of \jT\T ob- 
tained in the NMR shows a clear tendency against elec- 
tron doping. In the low doping region, which is near 
the magnetically ordered phase, one can see enhance- 
ment of \jT{F with decreasing temperature, while in 
the high doping region, there is no enhancement of 1/T\T 
and rather, it decreases with decreasing temperature. On 
the other hand, in the inelastic neutron scattering ex- 
periments, at least two common features have been ob- 
served: one is column-like response at the antifcrromag- 
netic wavenumber and another is gapped spin excitation 
spectra. 

Analysis of these results in the inelastic neutron 
scattering is often carried out with a Heiscnbcrg-type 
spin model, 25 ' 26 ) although this system shows itinerancy. 
Sometimes this system is regarded as a bad metal, but 



it is basically metallic in transport properties 27 -* (at 
most semiconducting). In addition, ARPES measure- 
ments show that the band structure coincides with those 
obtained in the first principle calculation with a mass 
renormalization factor of 2-4. 28 > 29 ) This means that, al- 
though electron correlation plays important role in the 
iron-pnictides, it is not so strong as in the high Tc 
cuprates, whose mother compounds are well described 
by a localized-spin Heisenberg model. Thus, in this pa- 
per we analyze the inelastic neutron measurements based 
on a purely itinerant model. We also discuss NMR results 
of \jT\T briefly to be more concrete. 

The model we use is the five-band Hubbard model 
proposed by Kuroki et al. 5 ^ which is down-folded from 
the first principle calculation. Five bands in this model 
mainly comes from five Fe-3d orbitals and hopping pa- 
rameters are kept up to fifth nearest neighbor. The dis- 
persion relation and the structure of the density of states 
are in good agreement with the first principle calculation. 
Although this model is two-dimensional, this simplifica- 
tion is not so harmful since the Fermi surfaces of this sys- 
tem are basically cylindrical. We must also note that this 
model contains one Fe atom per unit cell and Brillouin 
zone is doubled compared with the original one. (See up- 
per panel of Fig. 1.) Actually, this model is constructed 
from LaFeAsO system, but we believe that the proper- 
ties of FeAs layer is qualitatively same for other types 
of iron-pnictide. The interaction is limited to the onsite 
interaction and treated within RPA approximation. We 
tried several sets of interaction parameters keeping the 
constraints of U = U' + J+J 1 and J — J', but the results 
shown in the following are those with U = 1.1, U' = 0.9, 
J = 0.1 and J' = 0.1. Here, the unit of the energy is elec- 
tron volt (eV). We treat doping as a rigid band shift, i.e., 
we neglect possible structural change due to the change 
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of the charge balance between the conduction layers and 
blocking layers. Note that n — 6.0 corresponds to the 
stoichiometric sample in this model. We divide Brillouin 
zone into 64x64 meshes and 1024 Matsubara frequen- 
cies are used. Pade approximation is used to transform 
the results on Matsubara frequencies to those on real 
frequency. 

In order to analyze the inelastic neutron scattering 
experiments, we calculate spin structure factor S(q,uS), 
which is defined as Fourier transform of a spin correlation 
function, 



Ate 



l (S -S r (t)). (1) 



To be precise, we must include a form factor, often de- 
noted as F(q), in the above definition for the comparison 
with experimentally obtained scattering cross section. 
However, we ignore them since we believe that they do 
not change the following conclusion. In the spin-isotropic 
case, we calculate 



dTe^(StJr)S-) 



where is defined as 



C q+ka] Cka l' 



(2) 



(3) 



and a denotes the orbital indices. Then, S(q,u>) is ob- 
tained using the fluctuation-dissipation theorem as 

lmx + ~(q,oj) 



(4) 



1 - e-P" 

The main panel of Fig. 1 shows the overall view of 
S(q, lu) calculated at n = 6.0 and temperature T — 
0.02eV along the high symmetry lines of the Brillouin 
zone shown in Fig. 1 (a). Note that the energy range is 
beyond the experimental reach. We can see two promi- 
nent features in this figure: one is the peak at q ~ (tt, tt) 
and lu ~ leV and another is the peak at q ~ (tt,0) and 
lu ~ OeV. Considering the DOS (shown in the inset of 
the main panel of Fig. 1) and dispersion relations, we 
can see that the former corresponds to a transition be- 
tween two van Hove singularities above and bellow the 
Fermi energy. The latter peak comes from the nesting of 
the Fermi surface. Remember that the wave vector (tt, 0) 
in the present model corresponds to (tt, tt) in the original 
Brillouin zone (Fig. 1(b)). 

The low energy part of S(q, lu) is shown in Fig. 2 to 
compare with experimental results. We can see that the 
intensity distributions above (tt, 0) are almost vertical in 
this energy range, which is similar to the experimentally 
observed S(q,u>). Quite often, the column- like response 
at (tt, 0) is interpreted by localized spin model. However, 
the present results indicates that the localized spin model 
is not necessary to explain the spin dynamics of this sys- 
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Fig. 1. (Color online) Overall view of S(q, u) along the high sym- 
metry lines in the Brillouin zone at temperature T = 0.02eV. 
Inset shows DOS of the present model. The upper panels denote 
the Brillouin zone of the present model containing one Fe per 
unit cell (a), original Brillouin zone having two Fe per unit cell 
(b), and that in the magnetic phase (c). 



tern. In principle, S(q,Lu) in an itinerant model has low- 
lying excitations for < \q\ < 2k p. However, in the 
present model for iron-pnictides, there are disconnected 
small Fermi surfaces around (0,0) and around (tt, 0). In 
addition to this, since 2fcp is small, S(q,uj) has a rather 
narrow column-like response at (7r,0). Furthermore, the 
moderate strength of Coulomb interaction makes the re- 
sponse more sharper in RPA treatment. Actually, calcu- 
lated peak width is estimated to be about 0.125 of the 
distance between (0,0) and (tt,0) in the extended Bril- 
louin zone (In Fig. 2, S(q,Lu') with lu' = 15meV is plot- 
ted as inset). This width is in fairly good agreement with 
experimentally observed width. 24 ' Although the most re- 
ported data (not all 20 )) are lower than about 25meV at 
present, we make a plot of the same data up to 300meV 
in the upper panel of Fig. 2. We can see that the peak in- 
tensity vanishes with increasing frequency in this energy 
range and also that the V-shape peak structure becomes 
visible above O.leV. Here, in comparing the calculated 
result with the experiments, we must note that the band 
renormalization effect is neglected within our RPA treat- 
ment. Thus, the actual energy range will be lower than 
that indicated in Fig. 2, since the dispersion relation is 
squished by the band renormalization effect. Since the 
obtained peak intensity is so prominent at (it, Q), the 
result does not change if we take into account polycrys- 
talline nature. 

In order to obtain a clearer view, the w-dependence 
of S(Q,Ld) at Q = (tt, 0), is shown in the upper panel 
of Fig. 3. We can see that S(Q,lj) grows gradually with 
decreasing temperature. To see this behavior, the tem- 
perature dependence of the peak intensity, i.e. S(Q,lu') 
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Fig. 2. (Color online) Low energy close-up of S(q,u>) at temper- 
ature T = 0.02eV (upper panel: up to 300meV, lower panel: up 
to 40meV). Inset: S(q,a/) with ui' = 15meV. 



Fig. 3. (Color online) Frequency dependence of S(Q, u>) (upper 
panel) and Imx" 1 (Q,uj)/u) (lower panel) at Q = (n, 0) for sev- 
eral temperatures. Inset of upper panel: Temperature depen- 
dence of S(Q,u>') at Q = (tt,0) and w 1 = 15meV. Note that 
horizontal axis of inset is temperature (meV). 



at ui' = 15meV is plotted in the inset of Fig. 3. This 
result is similar to the experimentally observed temper- 
ature dependence of the peak intensity above the order- 
ing temperature. 23 ^ We also plot Imx H (uj,Q)/lu in the 
lower panel of Fig. 3 to compare it with the other ex- 
periment by McQueeney et al. 21 ^ This behavior is also 
consistent with the experimentally obtained result above 
ordering temperature, 21 ) considering the renormaliza- 
tion effect on energy scale. These results on tempera- 
ture and frequency dependence is usual in RPA theory 
of magnetism. Here, we make a technical remark. Shown 
frequency dependence is obtained by Pade approxima- 
tion, but we also calculate the same quantity in another 
method where analytic continuation is not needed. Two 
results are basically same, but Pade approximation gives 
a little smoother frequency dependence. 

Here, we briefly explain how above picture would 
change with doping. With electron doping, the peak posi- 
tion around u> — shifts to the incommensurate position, 
where the peak of the static susceptibility is found. 5 ' 
Namely, the main peak shifts in the direction from (7r, 0) 
to (tt, 7r). As expected, the scattering intensity of the low 
energy excitations quickly weakens by doping since the 
doping breaks the nesting condition and brings the sys- 
tem away from the magnetic instability. 

All the above results show that the itinerant model 
gives natural explanation for the neutron scattering ex- 
periments at least in the state without long-range or- 
der, although the sharp peak near (77, 0) has been often 
analyzed by local spin models. Here, we want to make 



some comments on the ordered state. First is about how 
S(q, uj) looks like in the ordered state. If we consider only 
Stoner type excitations as in the previous treatment, and 
if there appears an SDW gap in the ordered state, low 
energy electronic excitation is removed and S{q,Lo) in- 
tensity drops for to in the SDW gap. This scenario may 
explain the observed temperature dependence of peak 
intensity of S(q,uj) shown in Ishikado et al. 23 ^ However, 
in the ordered state, we must consider collective excita- 
tions, i.e., spin wave, as well. Then, it is not so simple to 
reveal the actual shape of S(q,ui) in the ordered phase. 
Second is about the order of the phase transition, i.e., 
second order or first order. In our treatment, or in the 
standard RPA approximation, transition becomes second 
order. However, in the actual compounds, lattice defor- 
mation is accompanied and this may bring the transi- 
tion to first order. Third, the consideration on the three 
dimensionality is required to discuss the ordered phase 
precisely. However, as shown by Matan et al. 24 ^ the spin 
dynamics is highly two dimensional above the ordering 
temperature, which justifies the present treatment. 

Finally, we calculate the NMR result of l/T x T on 
the same footing, i.e., itinerant model and RPA approx- 
imation. Figure 4 shows the temperature dependence of 
\jT\T obtained in the present treatment. Actually, Ikeda 
has calculated the normal state \ jT{T in the FLEX ap- 
proximation, in which the self energy correction is also 
taken into account, 30 ) and the arguments similar to our 
treatment are presented in Graser et al? x > Thus, the 
present calculation for \ jT\T is just for checking the pa- 
rameter choice. The obtained result in Fig. 4 basically 
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Fig. 4. (Color online) Temperature dependence of 1/TiT for sev- 
eral doping levels. 



captures the experimentally observed properties. 12, 19 - ) 
In the lower doping side, 1/TiT shows enhancement 
with decreasing temperature due to the Fermi surface 
nesting. On the other hand, in the higher doping side, 
\jT\T shows slight decrease, which comes from the dis- 
appearance of the Fermi surface nesting as well as the 
DOS structure of this model. However, this result is not 
enough to explain the experiments quantitatively. One 
is about the doping dependence. In experiments, the be- 
havior of \jT\T changes from increasing to decreasing 
more rapidly with doping. Another is about the behav- 
ior in the higher doping side. Observed decrease is more 
prominent than the present calculation, and shows one 
order drop from room temperature to superconducting 
transition temperature. Although the FLEX approxima- 
tion 30 ) gives a better result on this issue, the effects be- 
yond the rigid band shift may be important. 

In summary, it is shown that the itinerant model 
treated within RPA approximation can consistently ex- 
plain the inelastic neutron scattering results such as peak 
width or temperature dependence of the peak intensity. 
The NMR results can be also roughly reproduced on 
the same footing, but quantitatively, it is not satisfac- 
tory yet. Although our conclusion is that the itinerant 
model reproduces normal state spin dynamics well (at 
least qualitatively), our calculation does not exclude the 
possibility of the localized-spin model. Further works will 
be necessary for determining whether the strong corre- 
lation and/or Mott insulating behavior is important in 
this iron-pnictide superconductors. 
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